QUANTUM ERGODICITY OF BOUNDARY VALUES OF 
EIGENFUNCTIONS: A CONTROL THEORY APPROACH 



N. BURQ 

Abstract. Consider M, a bounded domain in W^, which is a Riemanian manifold with piece- 
wise smooth boundary and suppose that the bilhard associated to the geodesic flow reflecting 
on the boundary acording to the laws of geometric optics is ergodic. We prove that the bound- 
ary value of the eigenfunctions of the Laplace operator with reasonable boundary conditions 
are asymptotically equidistributed in the boundary, extending previous results by Gerard, Le- 
ichtnam T and Hassel, Zelditch 8: obtained under the additional assumption of the convexity 
of M. 

Resume. Soit M un domain borne de R'* qui est une variete riemanienne a coins. On suppose 
que le billard defini par le flot geodesique brise est ergodique. On demontre que les valeurs 
au bord des fonctions propres du Laplacien (avec des conditions aux limites raisonnables) sont 
asymptotiquement equidistribuees dans le bord. Ceci generalise des resultats anterieurs de P. 
Gerard, E. Leichtnam |7] et A. Hassel, S. Zelditch obtenus sous I'hypothese supplementaire 
de convexite geodesique du domaine. 



1. Introduction 

The purpose of this article is to show how ideas coming from control theory allow to deduce 
ergodicity of boundary values of eigenfunctions on manifolds with piecewise smooth boundaries 
and ergodic billiard flows from the known results on ergodicity of the eigenfunctions in the 
interior. The main result we obtain is the following: 

Theorem 1. Consider the eigenfunctions of the Laplace operator associated to the boundary 
condition B: 

— AoCo = A^e, 

(1.1) ^ ^ ^ ^ 

Bcj \qm = 

Denote by \dx'\ and da the natural measures on dM and T*dM ; and associate to ej its boundary 
value according to the table below 



B 


Be 


e' 


dfib 


Dirichlet 


e \dM 


\ ^dnC \qm 


7(f?)l|el,(.)<i^^ 


Neumann 


dne \dM 


e \dM 




Robin 


{dne - Ke) \3m 


e \dM 


7(f)) 1 ^ 
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where K e -^^{dM) is a non-negative self adjoint pseudodijjerential operator of order 1 (this 
assumption could be relaxed to the non negativness of its principal symbol). Then there exists a 
subset S of density 1 such that for any semi- classical pseudodifferential operator on the boundary 
(tangential to the boundary of the boundary) of order 0, a{x,hDx) we have 



The study of ergodicity of the eigenfunctions of the Laplace operator dates back to the works 
by Shnirelman ^Jj, Zelditch JJ| and Cohn de Verdiere [HI. In the case of manifold with bound- 
ary, the first study was initiated by Gerard and Leichtnam [Jj for C^'^ boundaries and then 
extended (with a much simpler proof) by Zelditch and Zworski to the case of manifolds 
with corners. 

Recently, the study of the ergodicity of the boundary values of the eigenfunctions has been 
performed in the restricted case of convex manifolds (with corners) by Hassel and Zelditch [H] . 
Theorem ^ is hence a generalization of these results in 1^ to the non convex case. Since most 
known ergodic billiards are non convex, our generalization is relevant (see the discussion in the 
introduction of [Hj). 

In fact it has to be noticed that for Dirichlet boundary conditions and the standard metric, 
the result above is not new (even though it does not appear explicitly in the literature for non 
convex domains): indeed it had been recognized by Gerard and Leichtnam [71 Remark 5.3] 
that for the Dirichlet problem, the questions of ergodicity of the eigenfunctions in the interior 
and ergodicity of their normal derivatives on the boundary are equivalent. The assumption of 
convexity of the domain which was made throughout the article played no role in this particular 
remark; and to obtain the ergodicity of the normal derivatives for a general ergodic billiard, it 
is consequently sufficient to link the posterior result in and this remark (and to treat some 
minor problems arising from the presence of the corners, see the end of Section 

On the other hand, in control theory, a lot of analysis has been done since the 80's to relate 
the behaviour of solutions of wave equations in the interior of a domain with the behaviour of 
their traces on the boundary (see 121 El ) • We are going to take advantage of this knowledge 
(eigenfunctions are -very- particular cases of solutions of wave equations) in order to show that 
the interior result in ^2] implies the boundary result. For general boundary conditions, the 
analysis becomes more involved than in the Dirichlet case; but the knowledge of (the eigen- 
functions are, for a density-1 sequence, equidistributed and consequently can not concentrate at 
the boundary) allows to avoid most of the technical difficulties (see |31) usually posed by bound- 
ary value problems; and we shall be able to perform quite the same argument as in 7 in the 
case of boundary conditions satisfying the strong Lopatinskii condition (non degenerate Robin 
conditions for example). In the case of a Neumann problem, this strong Lopatinskii condition is 
no more fulfilled at the gliding points (and in particular the traces of the eigenfunctions are not 
a priori bounded in L'^{dM)) and one has to show in addition that nothing happens at these 
points. In fact this part of the analysis, as well as the study of the behaviour near the singular 
part of the boundary, will be performed using some heat kernel methods from 8, Appendix]. 

In section 121 we give the assumptions and define the semiclassical measures associated to se- 
quences of eigenfunctions. We also give some elementary consequences of Zelditch and Zworski's 
result. We present an elementary and self contained (assuming the result in ^21) proof of The- 
orem n for Dirichlet conditions in section |31 In section 01 we prove the result for Robin and 



(1.2) 
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Neumann conditions for a restricted class of operators whose principal symbol vanishes near the 
singular set T, and the Glancing set Q. Finaly in sectional we explain how the analysis in [HI 
Appendix] allows to handle these singular and glancing sets. 

Remark 1.1. After this work was completed, S. Zelditch informed me that they had simul- 
taneously been able to supress the convexity assumption in their work with A. Hassel (see |S1 
Addendum] J. 

2. Assumptions and ergodicity of eigenfunctions 

The manifold (M, g) we consider is a compact C°° Riemann manifold with piecewise boundary: 
M is a compact subset of a C°° compact Riemann manifold, M and there exist r functions 
fi,...,fr£C°°{M) such that 

. |/ri(0)/0, 

• 2) dfi and dfj are independant on ff^{0) D fj^^{0), 

• 3) M = {xe M; fi{x) > OVi = 1, . . . , r} 

To include Bunimovitch's stadium example, condition 2) can be replaced by 

• 2') /j"^(0) n /j~^(0) is an embedded submanifold of M and M has Lipschitz boundary. 
A point xq G dM will be called regular if 

(2.1) tj{j;/^.(xo)=0} = l 

and singular otherwise. We will denote by S C dM the set of singular points. 

2.1. Ergodicity of eigenfunctions. With these assumptions, one can define the billiard flow 
(see p^) and under the assumption that this flow is ergodic, S. Zelditch and M. Zworski prove 
the following 

Theorem 2 (^U). There exist a set S Cf^ of density one: 

(2.2) hm h{Sn{l,2,...,N}) = l 

Af— >+oo iV 

such that for any A G "^phg pseudodifjerential operator of order 0, with kernel supported in the 
interior of M , one has 

(2.3) lim {Aej,ej)^2a.n = {dL,aQ{A)) 

j£o— >+oo ' 

where dL is the Liouville measure on S*M the sphere cotangent bundle. 

Remark 2.1. In jEj, this result is proved only for Dirichlet or Neumann boundary conditions. 
However, the proof relying on the easy (non-diffractive) part of Egorov theorem extends clearly 
to Robin boundary conditions. 
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2.2. Semi-classical measures. We introduce the following usual class of semi-classical symbols 
on M"': 

(2.4) S^{R^) = {a G C^{T*R^ x (0, 1]); \d^d^^a{x,^,h)\ < C,,/3,fc/i"^(0~l^l} 
To any symbol a G 5"^ we can associate an operator Op{a){x, hDx, h) by the relation 

(2.5) Op{a){x,hDx,h)u = j e'''-^a{x,hi,h)u{i)di. 

Denote by M the space M'^ endowed with a metric equal to g on M. Consider the sequence 
of eigenfunctions {ek)k&n normalized in L'^{M); and denote by their extension by outside 

M which is bounded in L?'{M). Using the standard definition of semi-classical measures in [7| 
section 2], |Uj (see also |3]) we can associate, modudo extraction of a subsequence, a measure 
on T*M such that for any semiclassical pseudodifferential operator Op(a)(a;, hDx) one has (with 
hk = ), 

(2.6) ^hrn^ (a(a;,/ifc-Dx)eA:,efc)^2(M) = 0)). 

Of course this measure is supported in T*M The semi-classical elliptic regularity implies 
that in the interior of M, the semi-classical measure fi is supported in the characteristic set: 

(2.7) supp(m) C Chaii-h^Ag - 1) = {{x, G T*R'^; \^\g(^x) = 1} 

The ergodicity property (|2.3() and 1)2. 7() imply that in the interior of M, /x satisfies 

(2.8) /xla;gA/ = ® %|g(^)=l- 

But due to the normalization ||ej||£2 = 1, we know that /U has at most total mass 1. Since it is 
also the case for dL ® <^|^|g(^)=i; we conclude 

(2.9) fi = ^IxeM = dL® l|5|g(^)=i 

and consequently for any semi-classical pseudodifferential operator on R*^, one has 



(2.10) . lim Op(a)(x,/i,Z),)e,-,e, =(dL®5|^l i,Op(a)(x,0>. 

Finally, we can extend slightly the class of symbols for which (|2.5|) holds. We work close to a 
point xq G dM \ S in a geodesic coordinate system (x„,x') such that locally M = {{xn-,x') G 
R*'+ X R'^-i} and 

(2.11) Ag = ^=dx„^/dJ)dx^+R{Xn,x\Dx') 

Vc(x) 

and the natural integration measures on M and dM are 



(2.12) \dx\ = c{x)dxndx' , \dx'\ = \J c{x)dx' 

In this coordinate system, the usual elliptic {£), hyperbolic (Ti) and glancing (Q) subsets of 
T*dM are 

£ =,{(x',r);i?(o,^',e' > 1} = {(^',e'); iriUx) > i} 

(2.13) g =,{(x',r);^(o,^',e' = 1} = {(^',e'); U'\\g(x) = 1} 
n =,{(x',0;m^',e' < 1} = {{^'.O; II^'IU.) < 1} 
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Lemma 2.2. Denote by S"^(Mx„ x T*R'^7^) the set of symbols in the x',^' variable depending 
smoothly of the parameter x„. For any oi, 02 G <S'''(Ma:„ x T*R'^7^) 

(2.14) lim {Op{ai){xn,x',hkDxi,hk)ek,ek)L^(M) = {fJ',ai{x,^'))- 

(2.15) lim (Op(a2)(x„,a;',/iA,.D^/)/iA,.D^„efc,eA..)^2(^) = (/^,a2(2;,^')Cn)- 

Indeed if the operators ai and 02 have kernels supported in M, this result is straightforward 
(using the elliptic regularity). To allow symbols supported up to the boundary, remark first that 
if 93 G C^(M) is equal to 1 close to 0, due to (|2.9() we have 

(2.16) lim 1™ yixn/e)ej\\l2^M) = lim{n,ip{-/ef) = 0. 
On the other hand 

= / {-h^jAg - l)ejip'^{x/e)-ej\dx\ 



(2.17) = (f (x / e)hj'V gejhjV gej\dx\ — ip {x/e)ejej\dx\ 

Jm Jm 



+ hj / Vg{(p'^{x/e))hjVgejej\dx\ + / h'^Ki 
Jm JdM 



(with = in case of Dirichlet or Neumann boundary conditions) which implies (using that 
the operator K is non-negative) 

(2.18) lim / ip'^ {x /e)hj\7 gCjhjV gej\dx\ < lim / ip'^{x/e)ej'ej\dx\ + lim Oe{hj) 
J^+'^Jm J^+'^Jm 

and consequently 

(2.19) lim lim / ip'^{x/e)hjVgejhjVgej\dx\ =0 

Now ()2.16l) . 1)2.19(1 (and the knowledge of the result if the operators are supported in the interior) 
imply Lemma 12.21 

3. The Dirichlet case 

3.1. A priori estimates. We start this section by recalling some well known estimates on the 
boundary values of the normal derivatives of the eigenfunctions (see for example in the more 
general context of strong Lopatinskii conditions |5| Theorem VII. 4. 4] or in our context ^ Lemma 
2.1]). 

Proposition 3.1. Consider (ej) the Dirichlet eigenfunctions of the Laplace operator associated 
to the eigenvalues Xj and normalized to 1 in L?'{M). Denote by e^ = Xj^dnCj \qm the normal 
derivatives on the boundary of the eigenfunctions. Then there exists C > such that 

(3-1) l|e5llL2{9M) < C 
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Proof. Consider X a vector field defined near M and such that on dM one has X ■ dn > O 0, 
and compute 



/ [-hlA - l,X]ekekdx = / -hlAXekek\dx\, 
Jm Jm 



(3.2) = K I Xck \dM -^\dx" 

^ ' 'dM on 



> c 







f 

IdM 


dn 



2 

\dx'\ 



On the other hand, the left hand term in (|3.1|) has the form 

(3.3) / hlP2ekej: 

Jm 



where P2 is a second order differential operator, and hence this quantity is clearly bounded by 

|2 

l//l(M) 



3.2. The relationship between the interior and the boundary. Consider a point xq € 
dM \ S. Working in a geodesic coordinate system close to xq, we can compute for symbols 
a = ai(x,^') + 02(2;, ^')Cn, using Green's formula 

(3.4) / -^[-hlAg-l,Ov{a){x,hkD,,)]ej-j\dx\ 
Jm ^i^k 

= I -7^{-hlAg - l)Op{a){x,hkDx)ejej\dx'\ 
Jm ^'ife 

= / ihkdx^Op{a){x,hkD:^)ejej\dx'\ 
JdM 

- / iOp{a){x,hkDx)ej \j:^=o hkdx„ej\dx'\ 
JdM 

= / OpOp(a2)(x„ = 0,x' ,hjDx>)hjdnej |^„=o hjdx„ej\dx'\ 
JdM 

Now consider the left hand side in (|3.4|) 
(3.5) 

/ —[-hlAg - l,Op{a){x,hkDx)]ejej\dx\ = / Op{q){xn,x' ,hD,j:„,hDx>)ejej\dx\ + 0{hj) 
Jm Jm 

where 

(3.6) q{xn, x',Cri, = {Cl + R{xn,x',^') - 1, a(x, ^)} 

is a second order polynomial in the ^„ variable. Taking into account the equation {—h^Ag — l)e = 
and (|2.11|) . we can eliminate the second order term in ^„ and obtain (using (|2.14|1 . (|2.15|l ) that 
the limit j £ S ^ +00 of the left hand side in (|3.4|) is equal to 

(3.7) {f,,q) 
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where q is the first order polynomial in the ^„ variable equal to q on the set {^'^+R{xn, x' , C') — ^ = 
0}. But the measure n is supported in this set and consequently we get 

(3.8) M = {f^,q) 
which implies that the left hand side in ()3.4I) tends to 

(3.9) -{fi,Hp{a)) = {Hp{fi),a) 

where p{x,(^) = (^.^ + R{x,S,') — 1. But since we know that the measure /i is equidistributed, we 
have far from the boundary Hp{^) = 0; and close to the boundary we can compute Hp{fj,) by 
the jump formula: 
(3.10) 

Hp(ii) = (2^n5x„(la;„>o)Ai = 2\/l - R5x^=o(giS^^^^q^'^\dx' \-2y/ 1 - Rdx„=o'Si6^^^_^q^(gi\dx'\ 
and we obtain using ()3.4() and ()3.9|) 

(3.11) lim / a2{xn = 0,x' ,hjDx,)hjdnej \x^=o hjdx„ej\dx'\ 

= [ a2Vl-RlR<i\dx'\\d('\ = (asVl - R, \dx\ 1b<i) 

JT'dM 

which implies Theorem ^ for operators whose kernel vanishes close to S. So far, working locally 
near a point which is not in S, we have only described the distribution of hjdnCi away from 
the singular set E. To get a complete description, we have to check that the norm can not 
accumulate at S. For this we are going to refine 1)3. 1|) . Close to a point xq & T, there exists /j^ 
and /jj (possibly non unique) such that fi^ixo) = fi2{xo) = 0. Take x ^ Co°(^) equal to 1 close 
to and consider 1)3. 1(1 with X replaced by 

P.12) X f !MOI±IZ«M) 



2 ,.l\fH(x)\ + \h,(x)\ 



and taking the limit i — > +oo (for fixed e > 0) we obtain 
(3.13) 



lim sup / X 

k^+oo JdM 



(/.,(x))2 + (/,,(x))2 



, dek 
on 



2 



<(^uiei^(.),x( ]x]) 



But taking into account that /i is equal to dL l|5|=i, we see easily that the right hand side 
in (|3.13|) is bounded by 0{e) because the symbol is bounded by 1/e but supported in a set 
of measure ©(e^). As a consequence, passing to the limit e — > 0, we obtain that the normal 
derivative can not accumulate at the singular set and consequently ()3.11l) holds for general 
semi-classical pseudodifferential operators on dM, tangential to E. 

4. Robin and Neumann conditions 

We consider in this section the eigenfunctions associated to the Robin or Neumann boundary 
conditions. 
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4.1. A priori estimates. By standard elliptic results (see for example |31 Appendix] in a 
slightly different context), we know that if ^ is a semi-classical pseudo-differential operator 
supported in the elliptic region of the boundary, <5, one has for any s S M, 

(4.1) \\Ael\\HsidM)^0{k^+oo) 

To deal with the hyperbolic region, we come back to (jH31). We obt ain using the boundary 
condition 

(4.2) / -^[hlAg-l,Op{a){x,hkD^)]ejej\dx\ 

0'2{xn = 0,x' , hjDxi){l — R + h^K'^){x' , hDx')ej'ej\dx'\ 



JdM 



Remark that in 1)4.21) due to the elliptic regularity in the boundary (|4.1() we can replace by 
^(x', hkD^>)K'^^{x', hkD^'Y with ^ equal to 1 for \(,'\g(a:) < 2 and ^ equal to for > 3. 

The second order pseudo differential operator h'^'i>{x' , Dx')K'^'^{x' , Dx')* can be viewed as a 
semiclassical operator of order and principal symbol {"^{x' ,(^')\'^ai{K)'^{x' ,(^'). 

The left hand side in (|4.2)) is a second order differential operator in the x„ variable. Using 
the relation ([2.11)) we see easily that in the left hand side of ()4.2j) we can eliminate the second 
order term in the normal derivative hkDx„ ■ Consequently the left thand side in ()4.2j) is bounded 
by C^(||efc||j;^2(j\,^) + hk\\'Vek\\L^(^M)) ^ 2C{A). On the other hand, for any q supported in the 
hyperbolic region, choosing 02 such that 

(4.3) -a2{xn = 0,x',eKl -R + ai{Kf){x,0 = \q\'' > 

we obtain, using sharp Carding inequality (and the fact that by the trace Theorem we know 
that ll^fc 11^^-1/2 < C which gives a control on the remainder terms). 



loc 



(4.4) WQekWrnu) <C- a2{xn = ^,x\hjDx>){l- R + h^ K^){x\hDx')ej U„=o eyjdx'l < C 

JdM 

Remark 4.1. Remark that the estimate above holds for any k, not only for a density-1 subse- 
quence 

4.2. The relationship between the interior and the boundary. We come back to ()4.2p . 
As in section 131 the left hand side in 1)4. 2|) tends (for a density-1 subsequence) to 

(4.5) -{f^,Hpia)) 

from which we deduce, taking the limit in the right hand side 

(4.6) lim / a2{xn = 0,x' ,hjDx',hj){l — R + h'^K'^){x' ,hDx')ej \x„=o'^\dx'\ 

= {VT^lR<i\dx'\\d^'\,a X {{1 - R) + \ai{K)f)) 
which is Theorem n for the restricted class of symbols which can be written under the form 

(4.7) b{x',a = axi{l-R) + \aiiK)\^) 
(which is true if b is equal to close to ^). 
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5. The singular and glancing regions 

So far we have obtained Theoremnfor symbols which vanish close to the singular and glancing 
regions. To deal with these regions, we are going to use the following result obtained from heat 
kernel considerations in |S1 Appendix] (not surprisingly since in this reference, the result was 
used for the same purpose): 

Proposition 5.1. For any e > there exists a function ips equal to 1 on Ti but with support 
sufficiently close to E such that 



(5.1) 



1 ^ 



Proposition 5.2. For any xq £ dM \ S and any e > there exists a symbol ^'^(x',^') equal to 
1 on G close to xq but with support sufficiently close to Q such that 



(5.2) 



1 ^ 



\{df,b,\^,{x',er)\<e 



Indeed the equality in Proposition 15.11 is [SI Lemma 7.1], while the inequality follows if the 
support of (f^ is taken close enough to S; and Proposition 15.21 follows from 8, (11.14)] and the 
Karamata Tauberian Theorem, as in the proof of |H1 Lemma 7.1]. 

To finish the proof of Theorem ^ for general operators, it is enough to prove that for any 
semi-classical pseudodifferential operator A on dM (tangential to S), one has 



(5.3) 



lim 



N 

E 



dM 



Op{a){x, XfD,, ^]^)e]e]da - {f,b, a{x, 0) 



By an easy argument, one deduces then for fixed A the existence of a density- 1 subsequence 
such that (|1.2|) holds, then by a diagonal argument used in jJOj) and the existence of a 
density- 1 subsequence such that (|1.2|) holds for any A follows. 

To prove (|5.3|) . we use first Proposition 15. II which allows to replace the operator A by an oper- 
ator A^ whose kernel is supported away from S, modulo an error of order ©^(e). Consequently 
it is enough to prove the result for operators whose kernel are supported away from S; and from 
now on we suppose that A satisfies this condition. 
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Write 



N 



N 



(5-^) 7vE / Op{a){x,X-'D,,X-')e';e^.da-{^,b,a{x,0) 

■ -1 "/ OAI 



^ ^ E I (Op(a)(x, XJ^D^) - Xj^D^,)Op{a){x, Xj^D,, Xj^)^e{x', Xj^D,,)* 



ey.da -{tib,a{x, 0(1 -\^ six', 



1 



N 



+ iV E / ^fD-,)Op{a){x, Xj'D,,Xj')^,{x', X-^^ D ..f e]e)da 

+ |(M6,a(x,e)|^.(x',Ol')| 
Since the operator 

(5.5) (Op(a)(x, AT^i?,, at!) - ^,(x', X-^^ D,,)Op{a){x, Xj^D^, Xj^)^,{x', Xj^D^,y 

has its kernel supported away from T, and its symbol equal to close to the glancing set Q, we 
obtain using (|I?T|) and that 



(5.6) 



dM 



Op{a){x,X-T^D^)e''je''jda - {fib,a{x,^)) 



is uniformly bounded with respect to j and tends to for a density-1 subsequence. Consequently 
the limit of the first term in (|5.4|) is equal to 0. On the other hand the contributions of the other 
terms in (jSHJ are, according to Proposition 15.21 of size 0{e). This concludes the proof of (|5.3j) . 
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